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SUMMARY ■ I 

The boundary conditions corresponding to the design problem when the blades being simulated by the 
bound vorticity distribution are presented. The 3D flow is analyzed by the two steps S2 — Si approach. In 
the first step, the number of blades is supposed to be infinite, the vortex distribution is transformed into an 
axisymmetric one, so that the flow field can be analyzed in a meridional plane. The thickness distribution of 
the blade producing the flow channel striction is taken into account by the modification of metric tensor in 
the continuity equation. Using the meridional stream function to define the flow field, the mass conservation is 
satisfied automatically. The governing equation is deduced from the relation between the azimuthal component 
of the vorticity and the meridional velocity. The value of the azimuthal component of the vorticity is provided 
by the hub to shroud equilibrium condition. This step leads to the determination of the axisymmetric stream 
sheets as well as the approximate camber surface of the blade. In the second step, the finite number of blades is 
taken into account, the inverse problem corresponding to the blade to blade flow confined in each stream sheet 
is analyzed. The momentum equation implies that the free vortex of the absolute velocity must be tangential 
to the stream sheet. The governing equation for the blade to blade flow stream function is deduced from this 
condition. At the beginning, the upper and the lower surfaces of the blades are created from the camber surface 
obtained from the first step with the assigned thickness distribution. The bound vorticity distribution and the 
penetrating flux conservation applied on the presumed blade surface constitute the boundary conditions of the 
inverse problem. The detection of this flux leads to the rectification of the geometry of the blades. 
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circulation 
upstream velocity 
pitch of the cascade 
inlet angle 
outlet angle 
potential function 
stream function 

bound vortex distribution function, or loading function 
tangential displacement 

camber line inclination angle with respect to the meridional plane 
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Superscripts 


angular velocity of the rotor 

pressure 

total pressure 

stagnation enthalpy or pt/p 
rothalpy or H + u)(V$) 
blade force 
dissipative force 
efficiency 

transformed coordinates system in Si approach 

leading edge 
trailing edge 
reference 

tangentiel component 
normal component 
nodal point indices 

upper side of the blade 
lower side of the blade 


1. INTRODUCTION. 

Most of the blading design procedures consider the velocity distribution on both sides of the blade as the 
initial data, the inverse problem becomes ill-posed and the designer loses the control of thickness distribution of 
the blade. To overcome this deficiency, this paper suggests an inverse method by representing the blades by a 
distribution of bound and free vortices which produce the desired swirl (V^r) variation. By introduction of the 
notion of associated elements on both sides of the blade in respect of the thickness distribution, and by imposing 
a conservative flux penetration through each pair of the associated elements when the geometry of the blade 
is not yet well defined, we obtain the well-posedness of the inverse problem. The iterative rectification of the 
camber surface in order to cancel the flow penetration leads to the final geometry of the blade. Treating first 
the 2D cascade design, §2 is devoted to show how to get the well posed inverse problem with the appropriate 
boundary conditions applied on the presumed blade contour, and the procedure leading to the rectification of 
the camber line related to the penetrating flux of the fluid determined on both sides of the blade. To treat the 
quasi 3D design, the S2 and SI approach as proposed by C.H. Wu [l] is adopted. The loading produced by the 
velocity difference between the two faces of the blade is directly related to the bound vorticity distribution that 
the blade has to generate. Assuming the number of blades infinite, the vortex distribution as well as the flow 
field become axisymmetric (S2 flow), §3 shows how the blade thickness distribution and the loading distribution 
can be taken into account in this scheme, and how to deduce the pressure distribution on the blades when their 
number is finite. An application to the case of the centrifugal impeller is presented. The loss scheme by the 
introduction of a plausible value of efficiency T} for each streamline as suggested by J.H. Horlock [2] is used. This 
approach opens up possibilities for the elaboration of a design which maintains the assigned value of the total 
pressure gain in each stage by modifying the (V^r) distribution in free space between blade rows. §4 is devoted 
to the blade to blade flow (SI) inverse problem, the boundary conditions for 2D inverse problem are transposed 
to this quasi-3D flow. The stream function is used to define the flow field and the finite volume method is used 
to solve the problem. Examples show the results concerning the design of centrifugal impeller. 


2. INVERSE PROBLEM FOR THE 2D CASCADE. 

Figure 1 shows the geometry of the blade characterized by its thickness distribution and the shape of 
its camber line. The arc elements taken respectively on the upper side and the lower side tangential to two 
inscribed circles centered on the camber line at x — dx/2 and x + dx/2 are called associated to the camber line 
element. The center of these associated elements are characterized by the abscissa x of the camber line element. 
Let V a represent the upstream velocity, h the pitch of the cascade, a and fi the inlet and outlet flow angles, the 
circulation T of the bound vortex generated by the blade is given by: 



167 


Third International Conference on Inverse Design Concepts and Opumization^n Engineering Sciences 
nrmF.S-Tin Editor. G.s Dnlikravich. Washin cipn D . C .. I77C - 


r = V 0 h(sin a — cos a tan 0) 

Th. bound vort.x distortion on ft. « « b. 

monotonic increasing function of x for the inverse problem. / ') ' A ^ b imposed near the trailing 

locul loading. Figur. 2 .how. ft. typwnl fonn of ft. ^,1 ft. ...» load.ng 

«*«' in o,d„ to obtain ft. ..,0 loading fording “, d , pled '„adi„g .dg. The how 

on both sides of the blade [3]. 


dip ..1+ TU^/ 


[4>]t = r/(x) or [j^ dl V- ~ T di dx 


( 2 . 1 ) 


A . the boundary condition is imposed on the 

be admitted., In order that the boundary condition does not produce extra 
each pair of associated elements is to be conservative, this implies: 


[?*dl] + _= 0 or Wt=0 


( 2 . 2 ) 


Bn 


The solution of the inverse problem determines the flux penetrating through the associated boundary elements, 
the camber line inclination correction 6d is given by. 

6 $ = 0.5[tan -1 (^jr) + + tan -1 (^) ] 

rxs^r ; S5sStt 

5. MERIDIONAL FLOW, S2 APPROACH. 

i first .tan the vortex distribution is transformed into an axisymmetric one by spreading it in the 

a, r tl^ 

stator is assumed to be oo, th. ... f the meridional velocity is 

Let f 1 f 2 = 6, and ( 3 represent the body fitted curvilinear coordinates (Fig. 4), the meriaio 

represented by: 0 = V»S» + V 3 *, = W't x + W 3 « 3 , the continuity equation becomes: 
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\fg simulating the elementary volume with striction in (3.1) is evaluated with § 22 - Using the stream function 
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y to represent the flow field by imposing: 
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the equation (3.1) is satisfied automatically. The governing equation for is obtained by writing V x V — n 3 e 2 , 
represents the azimuthal component of V x V^, it is deduced from the hub to shroud equ il ibrium 
condition. Let 
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( 3 . 3 ) 


In fact, there is a pressure gradient in the aximuthal direction in the flow space between blades, in the axial 
symmetric S2 flow where the number of blades is supposed infinite, this pressure gradient disappears and the 
volume force F b /p due to the blades has to be added in the momentum equation. The loss scheme [2] related 
to the plausible value of efficiency tj for each streamline of the stage is added, this scheme suggests that the 
dissipative force Ft/p i* related to the variation de V,r via q: 


f to - ,utor 

T ~ 1 (l-ri)jfrp\W VCMIVV rotor 


( 3 . 4 ) 


= o as well as F*t = 0 are imposed in the free space. Figure 5 shows the relation between the kinetic moment 
distribution in the blade row space and the circulation of the bound vortices produced by the blades. Let T+ 
denote the circulation generated by the blade in the section cut by an axis ym me trie stream surface \p etc, the 
kinetic moment (V*r) m , fc generated by the bound vortices located between the leading edge and the abscissa m 

can be represented by: 

(V # r) m> * = (Vir),..# + ^^/(m.V-) ( 3 5 ) 

Using (3.4) and adopting that dl/d£* or dH/d( 2 being equal to — ()2/P i® the dissipative scheme, the 
azimuthal component of the momentum equation leads to: 

(fi)a rtM^M ■ ,^Wi (3.6 

~~ [V ee se 1 


where W 3 = V t r + wr 3 and V 3 = V#r. The coordinates system C “ chosen so that the constant £ s lines are 
iteratively replaced by the streamlines. The component following t s of the momentum equation represents the 
hub to ahroud equilibrium condition, which gives: 
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Let n design the normal of the camber surface of the blade, we have: 


n = njF 1 -f njf 2 -f nsT* 


As W ± n in the rotor and V 1 n in the stator, we have: 


and Ft |( n, we have: 
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n 2 n 2 


(Ft)2 _ (Ft)i 


( 3 . 7 ) 


169 

Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 
nCTDES-nri. F-dnnr: G.S. D ulikravich. Washinnon D.C.. Qglflbr 2 3- 2 ^ . L° .2J - 
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In the free space, the component following T s of the momentum equation 

^ /-o 2 _ i r££ WO gWO 

free space y/gll ~yilges“ f2 g£3 


lead* directly to: 
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(3.8a) 
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The dot product of the momentum equation with V in the (tator and in the free .pace or with ft in the rotor 
lead* to the foilwing relation* which serve to update the nodal value* of E or I: 
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where 5( )/dm denote* the meridional streamwi*e tangential derivative. Writing V xV - ft 3 e 2 , we obtain 
the governing equation of ip: 


d / gu \ 5 / gas \ 

dp'py/gBt*** dp'py/gd 

_ d , 9n > 
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(3.10) 


For the inverse problem, the distribution of V t r is awigned, using (3.8), O i* updated iteratively. Let the 
camber surface of the blade be defined by « = ?{?,?) + etc, if the coordinate line, f 3 = ete are updated to 
the streamline* iteratively, ( 2 can be computed using the slip condition: 



(3.11) 


Figure 5 show, the geometry of the blading of a multistage turbopump obtained by solving the inverse problem. 
The CPU time on IBM 3090 in scalar mode is about 1 minute for the entire turbopnmp. The gnd used for ttae 
S2 computation is 300x16. Figure 6 shows the comparison of the centrifugal impeller* designed with 17 - 1 and 
t) < 1 having the same level of total pressure gain. 


Blade surface prea.ure evaluation. - Usually the S2 approach lead* to the determination of the mean 
velocity on both faces of the blade: 
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(3.12) 


Let AU denote the difference of the ab.olute velocities (V + - V) or the relative velocity [W - W ) «® the 
two face, of the blade, when the number of blade, is finite, this difference is related to the local density of 
bound vortex generated by the blade. In the S2 scheme, consider the blade section cut by * ( - «*« 
the flux of bound vortice. generated by the element 6? of the blade is determined by the J«* of fi through 

Z elementary surface (SS) 3 '> = 

that implies the circulation produced by AU is equal to the flux of the bound vertices we get the following 
r * Utl0D: , Am 2r co*ff (IV ,,+i/a.t, (3.13) 
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where denotes the local angle of the blade section with respect to the meridional plane. (3.12) and (3.13) are 
used to compute surface velocity on both faces of the blades, then the pressure distribution by the S2 approach 
can be deduced. (Fig. 9) 

4. BLADE TO BLADE FLOW, SI APPROACH. 

The blade to blade flow confined in each axisymmetric stream sheet is analyzed in order to define the final 
geometry for each section of the blade and to obtain the pressure distribution. At the beginning, the contour 
of the blade is created from the camber line obtained from the S2 step with the assigned thickness distribution. 
The conformal mapping (m, 8) => (x x ,x 2 ): 



(4.1) 


transforms the blade to blade flow confined in an axisymmetric stream sheet into a 2D cascade flow in the 
(x x ,x 2 ) plane. The body fitted coordinate system constituted by the equipotential lines f 1 = etc and the 
streamlines £ 2 = etc of a Active 2D flow around the cascade is created using the panel method [4]. In this 
system, the continuity equation becomes: 


jl [w^ vl) + = 0 


(4.2) 


where U x represent the contravariant components of the absolute velocity V for the stator and relative velocity 
W for the rotor and 

_ Dfc 1 ,! 2 ) , r .1 

where D(x 1 , x 2 )/D(£ 1 , £ 2 ) denotes the Jacobian, r represents the local thickness of the stream sheet. Introduc¬ 
ing the stream function tp with 
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(4.2) is satisfied. From the momentum equation, we can show that the free vortex of the absolute velocity 
shedding from the preceeding blade row must be tangential to the axisymetric stream sheet, the governing 
equation of the blade to blade flow stream function is deduced from this condition: for the relative flow around 
the blades of the rotor, we have: 


_[-*_( 

Ld£ lV 


9ii dtp. d_(9it_dtp_\\ = 
i/ + aci 


d^'pVgdt 1 ' ' 

dgiiW 1 , dg 12 W 2 


dC 


H" 




d log r 

» +2 ^t -ipt 


(4-4) 


Boundary conditions for the inverse problem: 

( Flux conservation: 


= o 


Bound vorticity assigned: [Widf 1 — c or^dO]^ = Tdf 


(4.5) 


The solution of the inverse problem leads to the determination of flux penetration on the blade contour, the 
camber line inclination correction 66 is given by: 

60 = 0.5[tan" 1 (^^-) + + tan_1 (^^)~] 


(4.6) 
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Figure 7 shows the network (£\ f 2 ) around a blade row for an impeller. Figure 8 shows the comparison 
of the camber lines of the impeller obtained from the S2 approach and rectified by the SI approach. Figure 
shows the pressure distributions obtained from the S2 and Si approaches. For the case of the turbopump the 
loading is optimised to avoid the cavitation. The results from the S2 and SI computations are similar, but not 
identical, the need of the SI computation to obtain the final geometry definition of the blades is confirmed, or 
one stream sheet, the CPU time on a IBM workstation RISC 6000/320 is about 40 minutes, or about 5 minutes 

on IBM 3090 in scalar mode. The grid used is 150x16. 


6 . CONCLUSION. 

The representation of the blades by the vortex distribution enables the formulation of the well-posed 
inverse problem, and which leads to design the blading of a turbomachine. The two steps S2 - SI quasi-3L) 
approach has been applied on different axial and radial geometries. Several kinds of loading function have 
been tried. The results show that the success of the blading design depends greatly on the meridional ( 9 r) 
distribution assignment associated with the loss distribution. To optimise the design in order to avoid the 
formation of the cavitation or the separation of the boundary layer in the design condition, when the loading 
is not too high, experiences show that an adequate modification of the bound vortices distribution function/ 
may effectively lead to prevent the surface pressure to be lower than the cavitation level or to maintain the 
adverse pressure gradient below the boundary layer separation criterion. The inverse problem procedure has 
been elaborated to calculate the turbomachines in incompressible range, the research works are planning o 
extend this method to make the transonic designs. 


ACKNOWLEDGEMENT 

Thi. work »* mpported by a ,..c«h contract from the SEP (So.Uti Europd.nn. d. Propnkion). The 
authors wish to thank them for their cooperation to this publication. 


REFERENCE 

1. Wu C.H., "A general theory of three dimensional flow in subsonic and supersonic turbomachines of axial, 
radial and mixed flow type. NACA TN 2604, 1952 

2. Horlock J.H., "On Entropy Production in Adiabatic Flow in Turbomachines.’’ Journal of Basic Engineering, 
Trans. ASME, Dec. 1971. 

3. Luu T.S., and Viney B., "The Turbomachine Blading Design Achieved by Solving the Inverse Problem. 
ASME Paper No.87-GT-215. 

4 Luu T.S., Coulmy G., and Viney B., 'Design problem of the profile or cascade of profiles and construction of 
the orthogonal networks using the Riemann surface for the multiform singularities." Computational Mechanic 

(1989) 4, 309-317. 



172 


Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 
flCIDES-Iin. Editor: G.S. Dulikravich. Washington D.C.. October 23-25. 1991._ 





Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 


rtmijsrauma 


UWl 


[IMflTFiSml 


Washington 


■r 23-25. 199 











174 


Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 
nODES-IID. Editor: G.S. Dulikravich. Washington D.C.. October 23-25. 1991. _ 



*7 = 


1 





t ) = 0.7 

Fig. 6. Centrifugal impellers designed with t] = 1 and r] < 1 
having the same level of total pressure gain. 
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